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ABSTRACT 



Aims. We re-examine the recombination/collisional emission line (RL/CEL) nebular abundance discrepancy problem in the light of 
recent high-quality abundance determinations in young stars in the Orion star-forming region. 

Methods. We re-evaluate the CEL and RL abundances of several elements in the Orion nebula and estimate the associated uncer- 
tainties, taking into account the uncertainties in the ionization correction factors for unseen ions. We estimate the amount of oxygen 
trapped in dust grains for several scenarios of dust formation. We compare the resulting gas+dust nebular abundances with the stellar 
abundances of a sample of 13 B-type stars from the Orion star-forming region (OriOBl), analyzed in Papers I and III of this series. 
Results. We find that the oxygen nebular abundance based on recombination lines agrees much better with the stellar abundances than 
the one derived from the collisionally excited lines. This result calls for further investigation. If the CEL/RL abundance discrepancy 
were caused by temperature fluctuations in the nebula, as argued by some authors, the same kind of discrepancy should be seen for the 
other elements, such as C, N and Ne, which is not what we find in the present study. Another problem is that with the RL abundances, 
the energy balance of the Orion nebula is not well understood. We make some suggestions concerning the next steps to undertake to 
solve this problem. 
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1. Introduction 

The knowledge of the chemical composition of stars and of the 
interstellar matter is the basis of our understanding the chemi- 
cal evolution of galaxies. In particular, elements such as carbon, 
nitrogen, and oxygen provide vital clues for the stellar evolution 
and nucleosynthesis as well as on the hi story of galaxies (see e.g . 
proceedings of CNO in the Universe bv lCharbonnel et aUl2003l) . 
H ii regions have for a long time been privileged sites to mea- 
sure chemical abundances in galaxies because they were easily 
observable at large distances and because the techniques to an- 
alyze their spectra in terms of abundances are simple. However, 
the last two decades have shown the important potential of stars, 
now that computers allow the calculation of very elaborate stel- 
lar atmosphere models and large telescopes allow one to obtain 
suitable spectra of individual stars not only in the Milky Way, 
but also in other galaxies. Indeed, stars are superior to nebulae 
as abundance indicators in that they are not affected by deple- 
tion into dust grains. Another problem with H n regions is that 
depending on the methods used to derive abundances, the results 
are different. When using collisional emission lines (CELs), the 
derived abundances are system atically smaller than when using 
recombination lines (RLs) (see iGarcla-RoTa s & Este banl I2007I 
and references therein). What is the reason for these discrep- 
ancies? And which one of the two abundances (if any) should 
one believe? After sever al decades of stu dies on this, the de- 
bate is still not closed. iPeimberJ ( 1 1967b has argued for the 
presence of temperature fluctuations (of still unknown cause) 
that bias the nebular abundances derived by CELs downwards. 
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A recent reassessment of t h is point of view can be found in 
iPeimbert & Peimbertl d2009l) . IStasihska et al.l d2007l) have exam- 
ined the hypothesis of inhomogeneous abundances in the inter- 
stellar medium (ISM) and concluded that if this were the case, it 
is the RL abund ances which would be severely biased upwards. 
lErcolanol (120091) showed that CEL-RL abundance discrepancies 
can be obtained in a chemically homogeneous medium if it con- 
tains X-ray-ionized dense clumps. However, she finds that the 
required X-ray flux to accound for the observed abundance dis- 
crepancy in the Orion nebula far exceeds the Orion's known stel- 
lar and diffuse X-ray budget. 

One way to attack this problem of abundance bias is to con- 
front abundances derived in Hn regions with those derived in 
B-type stars in the same environment, since both type of ob- 
jects are expected to share the same chemical composition, re- 
flecting the present-day abundance pattern of the galactic re- 
gion where they are located. This has been done, for example 
in the galaxy NGC300, where oxygen abundances have been 
derived for blue supergiant s as well as for giant H n regions by 
iBresolin et al.l (120091 120101) . The conclusion of that study was 
that stars and H n regions (CELs) give overall the same abun- 
dances. However, in view of the existence of important abun- 
dance s catter in galaxies (demonstra ted for the spiral galaxy 
M33 by iRosolowskv & SimonL 120081) . more accurate tests are 
needed. In particular, one should compare the abundances of B- 
type stars with those of the associated Hw region. 

In the Milky Way, iPrzvbilla. Nieva & Butleil (120081) com- 
pared the abundances of a small but representative sample of 
B-type stars in the solar v icinity with the abun dances in the 
Orion nebula published by lEsteban et alJ (120041) and with the 
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abundances in the neutral ISM from absorption line measure- 
ments. They concluded that taking into account depletion in 
dust, all sets of abundance s are r emarkably similar. The abun- 
dances from lEsteban et al.l d2004l) were those derived from re- 
combination lines, which would argue in favor of RLs giving the 
right nebular abundances. An update of the study of Przybilla 
et al. was however warranted, with a more critical analysis of 
CEL and RL abundances in the Orion nebula and a more com- 
plete stellar data base centered in the Orion star-forming region 
(Ori OBI), to better adress the RL/CEL abundance discrepancy. 

In ISimon-Dfazl d2010L Paper I), the first of this series 
of papers, we us ed the stellar atmosphere code FASTWIND 
dPuls et al.Ll2005l) to perform a thorough self-consistent spectro- 
scopic analysis with a high-quality set of spectra of 13 early B- 
type stars from the various subgroups comprising the Orion OB 1 
association. In this study it was shown that the dispersion of 
O and Si abundances between stars in the various subgroups 
found in previous analyses (e.g. lCunha & Lambertl 1 199211 19941) 
was a spurious result and the consequence of a bad character- 
ization of the abundance errors propagated from the uncertain- 
ties in the stellar parameter determination. This result is con- 
firmed in Nieva & Simon-Diaz (in prep., Paper III) by means 
of an independ ent spectroscopic analysis with the Atlas9 code 
dKuruczl [19931 and recent versions of DE TAIL and SURFACE 
iGiddingslligsTtlButler & Giddingslll985h . In paper III we also 
extend the elements analyzed to C, N, Mg, Ne, and Fe. 

In this paper, we focus on the stellar vs. nebular abundance 
comparison taking into account the effect of depletion onto dust 
grain in the nebula. The paper is structured as follows: In Sect. 2 
we present a brief summary of the abundances in B-type stars in 
the Orion OBI association from our recent studies. In Sect. 3 we 
present a critical analysis of CEL and RL abundances of various 
elements in the Orion nebula. In Sect. 4 we analyze the depletion 
of the elements in the dust phase that is occuring in the Orion 
nebula. In Sect. 5 we draw inferences from the comparison of 
stellar and nebular abundances in Orion. The main conclusions 
of this work are summarized in Sect. 6. 



2. Abundances in B-type stars 

Table Q] summarizes the stellar abundances derived in Papers I 
and III. Two types of uncertainties are given along with the mean 
values. The first one refers to the dispersion of abundances ob- 
tained for the 13 analyzed stars. The second one (in parenthe- 
ses) indicates the uncertainties intrinsic to the method of abun- 
dance determination in a given star (i.e. accounting for the effect 
of stellar parameters, microturbulence, and line-by-line abun- 
dance dispersion^]). The derived abundances agree with those 
computed by iPrzvbilla. Nieva & Butler] (120081) for a small sam- 
ple of B-type stars in the solar vicinity, confirming the homo- 
geneity of present-day abundances in the solar neighborhood as 
derived by young massive stellar objects already suggested by 
these authors. While the intrinsic uncertainties are on the order 
of 0. 1 dex, the abundances obtained in the various stars of our 
sample are remarkably coherent and lead to an unprecedented 
accuracy of the overall stellar abundance pattern in Orion. 



Table 1 Summary of abundances in B-type stars in Ori OB 1 . Data from 
ISimon-Dfazl d2010t Paper I) and Nieva & Simon-Diaz (in prep., Paper 
III). Mean value and standard deviation of abundances obtained for the 
sample of stars. Values in parentheses indicates the intrinsic uncertain- 
ties of the analysis of individual stars. 





e(X) 




6(X) 


c 


8.35 + 0.03 (0.09) 


Si 


7.5 1 1 ' +0.03 (0.08) 


N 


7.82 ± 0.07 (0.09) 


Mg 


7.57 ± 0.06 (0.03) 


O 


8.74^0.04(0.10) 


Fe 


7.50 + 0.04(0.10) 


Ne 


8.09 + 0.05 (0.09) 







' Mean value from both studies. 



3. Gas phase abundances in the Orion nebula 

The gas phase abundances in an ionized nebula can be read- 
ily obtained from observed intensity ratios of its emission lines 
through methods that are w ell-known and described e.g. in 
lOsterbrock & Ferlandl(l2006l) . After correcting the observed line 
fluxes for interstellar extinction, determining the temperature 
and density of the emitting zones with appropriate line ratios, 
the ionic abundances are obtained directly from emission line 
ratios. Elemental abundances are then obtained by summing the 
abundances of the observed ions and accounting for unseen ionic 
species with ionization correction factors (icfs). The chemical 
composition determined in this way depends on the quality of 
the observational data, the adopted dereddening procedure, the 
temperature, density indicators, and the chosen icfs. The proce- 
dures are quite standard, but may deviate in some detail from 
one author to another. Of course, the results depend on the set of 
atomic data that are used. 



3.1. O 

O from CELs Oxygen is the element for wich the abundances 
can be determined most reliably in Hn regions, because both 
ions that are expected to be present, + and 2+ , emit strong 
lines in the optical and, in addition, the temperature and density 
of the emitting zones can be easily obtained from appropriate 
line ratios. 

There have recently been a number of optical observations of 
the Orion nebula with high resolution and very good signal-to- 
noise ratio. To obtain a set of abundances in the Orion nebula that 
gives a fair account of uncertainties and possible variations from 
one point to another, we rederived in a consistent way the CEL 
abundances of oxygen in different positions of t he Orion neb- 
ula fro m the reddening-cor r ected intensi ties of by l Esteban et al.l 
(I1998L E981. lEsteban et all (12004 E04), iBlagravel (120061 B061 
iMesa-Delgado et alj d2009l M09). We included the two slit po- 
sitions provided by E98, B06, and M09. In the case of E98, 
both spectra correspond to the nebula, while for B06 and M09, 
one of them is located in the shock front of one of the Herbig- 
Haro objects of the Orion nebula. For the abundance analysis, 
we assumed that the electron density throughout is the one given 
by the [On] 3726/3729 line raticQ We have considered that H/J 



1 For Mg, the corresponding value does not account for line-by-line 
abundance dispersion, because only one line is available for the abun- 
dance analysis. 



2 Densities given by [CI m] or [Ar iv] doublets would, in principle, be 
more appropriate for the high-excitation zone, but they are much more 
uncertain in the Orion nebula. Fortunately the abundance analysis is 
not strongly affected by the density in the expected range of densities. 
Densities derived from the [Sn]67 17/6731 line ratio are slightly lower 
than those derived from [On]3726/3729. 
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Table 2 Newly computed densities, temperatures, and ionic abundances from optical spectra of the Orion nebula. Quartiles are indicated in small 
fonts. The references for the corresponding observations are indicated in the text. 



Observ. 


A'e(OII) 

[iu cm j 


r e (om) 

[10 3 K] 


r e (NH) 

[10 3 K] 


e(0 + ) 


e(O z+ ) 


6(0) 


E04 


5.10 

" 4.72 


8.31 


8.33 
S.2S 


10.3 


10.5 
10.1 


7.76 


7.80 
7.72 


8.44 


8.44 
8.43 


8.52; 


3.53 
i.51 


E98 (1) 


3.68 1| 


8.21 


8.63 
7.71 


10.1 


10.7 
8.97 


7.70 


7.98 
7.57 


8.41 


8.56 
8.26 


8.49 


3.67 
3.37 


E98 (2) 


6.08 ™ 


8.15 


8.37 
7.69 


10.9 


12.1) 
9.01 


7.79 


8.22 
7.59 


8.46 


8.62 
8.35 


8.54 


3.82 
S.47 


M09 (N) 


2 17 2 37 


8.15 


8.30 
7.95 


9.65 


9.94 

9.3S 


7.97 


8.02 
7.92 


8.37 


S.4I 
8.33 


8.51 


3.55 
3.48 


M09 (HH) 


22-00 £° 


8.73 


8.85 
8.49 


8.91 


9.40 
7.99 


8.42 


8.75 
8.23 


8.11 


8.18 
8.08 


8.59 


3.85 
3.46 


B06 (N) 


4.05^ 


8.48 


8.51 
8.46 


10.0 


10.0 
9.96 


7.87 


7.87 

7.86 


8.39 


8.40 
8.38 


8.51 


3.51 
350 


B06 (HH) 


Q 28.6 

4.36 


8.41 


8.47 
8.11 


9.80 


10.5 
7.49 


7.51 


8.27 
7.27 


8.70 


8.79 
8.68 


8.73 


3.91 
3.70 



and [Oiii]5007 are emitted in a zone whose temperature is ob- 
tained from the [O in]4363/5007 line ratio and that lines from 
[On]3727 and [Nn]6584 are emitted in a zone whose temper- 
ature is obtained from the [Nn]5755/6584 line ratio. The refer- 
ences for the atomic data used in the calculations are listed in 
Table [7] in Appendix A. The intensities of the forbidden lines 
were obtained from a 5-level model atom. 

The resulting electron densities, temperatures, and abun- 
dances of + , 2+ , and O in units of e(X)= 12 + log X/H are 
listed in Table [2] The uncertainties (listed as quartiles in the ta- 
ble) were obtained by Monte-Carlo simulations, taking into ac- 
count observational errors in the dereddened line intensities as 
listed in the original papers^. 

We find that CEL results are remarkably consistent with each 
other once the same atomic data and procedures are applied to all 
the spectra, when disregarding the two Herbig-Haro objects (la- 
belled HH in the Table). This is because of the exceptional qual- 
ity of the data reported in the original papers. As for our B stars, 
the dispersion of the derived CEL oxygen abundances is smaller 
than the formal errors: excluding the two Herbig-Haro object, 
the average value of e(O) is 8.51 and the dispersion is 0.02 dex. 
For convenience, below we will exclusively use the abundances 
derived from the spectra by E04. 

O from RLs The method using recombination lines to deter- 
mine ionic abundances has the advantage of being almost inde- 
pendent of errors in the electron temperature. On the other hand, 
it relies on very weak lines and requires a certain know-how to 
choose the best strategy for the different observed multiplets. We 
therefore adopted the ionic abundances derived in the original 
paper. Note that the abundances derived by E04 from individ- 
ual recombination lines of the same ion differ much more than 
expected from observational errors only. The final adopted abun- 
dances are a weighted mean of the sum-values of the most reli- 
able multiplets. The error attributed to the final RL abundances 
by E04, which we adopt here as well, is also a result of this 
weighting process. 

lEsteban et al.1 (120041) derive RL abundances for both + and 
2+ . But they note that the RL + abundance, which is based on 
a single faint line in a spectral region much affected by sky emis- 
sion, is not very reliable. We will thus use only their abundance 
of 2+ , and from it derive the RL oxygen abundance using the 



3 We have taken into account the fact that electron densities are ob- 
tained from adjacent lines whose measured intensity ratios do not de- 
pend on calibration errors, but only on line-fitting and signal-to-noise 
ratios. 



O /0 + obtained from forbidden lines. The result is presented in 
Table[3] 

3.2. N, Ne, S, Ar 

N, Ne, S, Ar from CELs We use the same procedure as for oxy- 
gen to derive the CEL ionic abundances of N + , Ne 2+ , Ar 2+ , S + , 
and S 2+ from the E04 spectrum, using the atomic data listed in 
Table [7] Like for oxygen, the uncertainties in the ionic abun- 
dances are determined from Monte-Carlo simulations of the ob- 
servational errors. The resulting ionic abundances are displayed 
in Table [4] To obtain the corresponding elemental abundances, 
we need to use ionization correction factors. Previous studies 
(E98, E04) used icfs from a variety of sources that are not nec- 
essarily consistent and sometimes outdated. To obtain the most 
reliable set of icfs for our purpose, we constructed a tailored 
photo ionization model of the Orion nebula with the Cloudy 
code dFerland et all 1 19981) . The way we proceeded to obtain the 
model is described in detail in Appendix [B] From this model 
we derive icfs for N, Ne, S and AiQ Those and the resulting ele- 
mental abundances are listed in Table|4] There are non negligible 
differenc es between our total abundances and the ones found by 
E04 (and lGarcfa-Rojas & EstebanLl2007t) . They stem essentially 
from the different icfs adopted in each case. Our values should be 
more reliable, because they were obtained from a tailored model. 
Note that some icfs are large, especially that of N + , which we 
find to be on the order of 10 (E04 attributed a value of about 
6). In this case, even with a tailored photoionization model, we 
cannot be entirely confident in the result we find. Concerning 
neon, the icf derived from our model implies that Ne/O is not 
equal to Ne ++ /0 ++ , which is the classical recipe to derive the 
neon abundance. Using th e spectroscopic data for a bout 80 zones 
in the Orion nebula from lO Dell & Harrisl d2010l) . we find that 
Ne ++ /0 ++ shows no trend with ++ /H over a range of one decade 
in ++ /H. This would rather support the classical icf for Ne over 
the one we derived from our model. Indeed, model atmospheres 
become more uncertain at energies above 40 eV, and the neb- 
ular constraint provided by our photoionization m odel is given 
by the very weak [Ariv] lines. On the other hand. iRubin et al.l 
(1201 Oh used infrared line intensities measured by the Spitzer 
Space Telescope to derive 12 + log (Ne/H) = 8.00 + 0.03 with- 
out the need of any icf because both [Ne n] and [Ne in] lines are 
observed. This value is very close to the one we derive with our 



4 Ar and S abundances have not been determined so far in our sample 
of stars, so they will not be used in this study. However, the stellar values 
should be available in the near future. 
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Table 3 RL abundances of O, C, and Ne from the E04 optical spectrum of the Orion nebula. 





Orl 


Crl 


Ne RL 


ion abundance 


e(O 2+ )=8.57±0.01 


e(C 2+ )=8.34±0.02 


e (Ne 2+ )=7.95±0.09 


icf 


log icf (O 2+ )=0.08±0.02 


log icf (C 2+ )=0.03±0.01 


log icf (Ne 2+ )=0.29±0.02 


element abundance 


6(O)=8.65±0.03 


6(C)=8.37±0.03 


e(Ne)=8.24±0.11 



Table 4 Newly computed CEL abundances of N, Ne, Ar, S from the E04 optical spectrum of the Orion nebula. 





NcEL 


Ne C EL 


Ar CE L 


Scel 


ion abundance 


e(N+)=6.92±0.05 


6(Ne 2+ )=7.76±0.01 


e(Ar 2+ )=6.37+0.02 


e(S 2+ )=6.80+0.05 


icf 


logicf(N + )=1.00±0.04 


log icf (Ne 2+ )=0.29±0.02 


log icf (Ar 2+ )=0.02±0.01 


log icf (S 2+ )=0.07±0.01 


element abundance 


e(N)=7.92±0.09 


6(Ne)=8.05±0.03 


e(Ar)=6.39±0.03 


e(S)=6.87+0.06 



model icf. In view of this dilemma, we give in Table[5]the values 
of Ne derived both with our model icf and with the classical icf. 



Ne from RLs Of the elements N, Ne, S and Ar, reliable abun- 
dances from recombination lines exist only for Ne, using the 
Ne 2+ ion. The results are listed in Table [3] again using both the 
model icf and the classical one. 



3.3. C 

For carbon there are no CELs in the optical wavelength range, 
and one has to resort to ultraviolet observations. There are sev- 
eral difficulties with this: the extinction correction is important 
and may be inaccurate. Collisionally excited ultraviolet lines de- 
pend very strongly on the electron temperature, so using temper- 
atures derived from [N n] or [Om] lines may lead to substantial 
errors. iRubin et"al] (1 19931) performed a detailed photoionization 
modeling to infer the carbon abundance by direct comparison of 
model results with observational data in several locations of the 
Orion nebula. They point out however that their model was not 
able to explain the Cm] 1907,09 and the Cn] 2324-29 lines at 
the same time. As a result, their carbon abundance is not very 
accurate. We will use the value e(C) = 8.4 and assign to it an 
uncertainty of +0.25 dex 

In the optical, there are several RLs that allow one to de- 
termine the abundance of C 2+ . We use the value determined by 
E04 and list it in Table [3] To determine the total carbon abun- 
dance from RLs, we use the icf derived from our photoioniza- 
tion model, as done for the other elements. The result is listed in 
Table [3] 



3.4. Fe 

There are many forbidden lines of Fe + , Fe 2+ and Fe 3+ in the 
optical spectra of the Orion nebula. To make the best use of 
them, multilevel atoms of over 15 levels have to be considered. 
iRodrfguez & Rubinl d2005l) performed a detailed analysis of the 
Fe spectrum in the Orion nebula using the spectrum of E04, and 
compared the Fe abundance obtained by summing the observed 
abundances of Fe + , Fe 2+ and Fe 3+ with the abundance derived 
using an icf implied by their photoionization model analysis. 
They found an inconsistency of about a factor 4, implying that 
one or several of the atomic data involved in the analysis are 
wrong. 



3.5. Mg and Si 

Collisionally excited lines of Si and Mg ions are not present 
in the optical range. Observations are available in the UV: 
Si ml 1883,93 and Mg n2768. The former have been analyzed 
bv lRubinetal] d!993l) . together with Cm] 1907,09. These au- 
thors find a Si/C ratio of 0.016 with an accuracy of about 30%. 
The absolute abundance of Si is plagued however with the uncer- 
tainty in the C abundance as deduced from UV lines. The abun- 
dance of Mg was inf erred from a photoionization modeling by 
iBaldwin et all d 19911) . but its value is not accurate, because the 
observed intensity of Mg u 2768 is very uncertain. In addition, 
the flux of Mg ii 2768 received at Earth is affected by interstellar 
absorption. 

3.6. Summary of gas phase abundances 

In Table[5]we summarize our nebular abundance analysis by giv- 
ing our adopted values of the gas phase C, N, O, Ne, S, Ar, Mg, 
Si, and Fe abundances. The error bars now include the uncer- 
tainties from the Monte-Carlo analysis and the uncertainties in 
the icfs, all added in quadrature. The element with the best deter- 
mination is oxygen. For nitrogen, the large icf introduces an un- 
certainty difficult to evaluate, but possibly larger than indicated 
in the table. For neon, the icf derived from our photoionization 
model leads to an abundance larger by 0.2 than the one obtained 
using a classical icf. This issue needs to be investigated further. 
The abundance of carbon from CELs is quite uncertain, as are 
the abundances of the refractory elements Fe, Si and Mg, each 
one for different reasons. 

Figure Q] compares the abundances of C, N, O, Ne, Mg, Si, 
and Fe in the Orion nebula obtained using CELs (green boxes) 
and RLs (red boxes) with the abundances in the B stars of 
Orion OBI (cyan boxes). In each case, the height of the boxes 
represents the uncertainties, indicated in Table Q] (standard devi- 
ation from the 13 stars) and Table [5] (Orion nebula). 

4. Depletion of the elements in the dust phase 

Before comparing the nebular abundances with the stellar ones, 
one has to consider that some of the material in the nebula may 
be in the form of dust grains, whose total mass and chemical 
composition cannot be directly accessed. That the abundances 
of the refractory elements Mg, Si, and Fe in the gas phase of the 
Orion nebula (Table |5]l are lower than those found in the atmo- 
spheres of our sample of B stars (Table[T|) by at least one order of 
magnitude clearly indicates that these elements are mostly found 
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Table 5 Summary of gas-phase abundances in the Orion nebula. 



A 


6(A) 


data 


source 


method 


major uncertainties 


C 


8.40+0.25 


CEL: Cm] 1907,09 


Rubin etal. (1993) 


model 


dereddening, CII not fitted 




8.37+0.03 


RL: Cn 


C ++ : E04 


icf (own model) 




N 


7.92+0.09 


CEL: [Nn]6584 


data: E04, N + : this paper 


icf (own model) 


icf 


O 


8.52 ±0.01 


CEL: [Om]5007, [On]3727 


data: E04, + & ++ : this paper 


no icf needed 






8.65+0.03 


RL: Om 


++ : E04 


icf (own model) 




Nc 


8.05+0.03 


CEL: [Nem]3869 


data: E04, Ne . this paper 


icf (own model) 






7 84+0 03 


CEL: [Nem]3869 


data: E04, Ne ++ '. this paper 


icf (classical) 






8.25±0.35 


RL: Nen 


Ne ++ : E04" 


icf (own model) 






8.03±0.26 


RL: Nen 


Ne ++ : E04" 


icf (classical) 




S 


6.87+0.06 


CEL: [Sm]9069 


data: E04, S ++ : this paper 


icf (own model) 




Ar 


6.39±0.03 


CEL: [Arm]7135 


data: E04, S ++ : this paper 


icf (own model) 




Mg 


6.50:: 


CEL: Mg ii 2798 


Baldwin etal. (1991) 


model 


intensity, interstellar absorption 


Si 


6.50+0.25 


CEL: Sim] 1883,92 


Rubin etal. (1993) 


model 


Si/C more accurate than Si/H 


Fe 


6.0+0.3 


CEL: [Fen], [Fern], [Few] 


Rodnsuez & Rubin (2005) 


model 


sum of observed ions gives 



different result 



"The uncertainty for Ne ++ was erroneaously given as 0.09 dex in E04, instead of 0.2 dex (J. Garcfa-Rojas, private communication). 



O 



_] 

. I_J _l - 
u 01 


_J 

_ LJ 1 _ 

o 01 


_i 

i_i_i i 

CJ 01 


_i 

. i_j _i _ 

O 01 


_l 

_ LJ 1 , 

O 01 


_i 

i_i_i i 

U 01 


_] 

_ U _l - 
O 01 


1- 






































1 
















1 


c 


N 





Ne 


Mg 


Si 


Fe 



Table 6 Dust composition possibilities considered in this study. Values 
in parentheses correspond to the fraction of molecules from each type. 



Fig. 1 Comparison of C, N, O, Ne, Mg, Si, and Fe abundances 
obtained from B-type stars (cyan), CEL (green), and RL (red). 
For Ne there are two green boxes, corresponfing to different icfs 
(see text). The height of the boxes represents the uncertainties, 
indicated in Tables Q]and[5] Note that for Mg there is no evalua- 
tion of the uncertainty. See text for more detail. 



in dust grains. We find, by comparing the numbers in those two 
tables, that ~ 90% of Mg, ~ 80-90% of Si, and ~ 93-98% of Fe 
are trapped in dust grains. If we knew the chemical composition 
of dust, we could estimate the degree of depletion of the other 
elements. 

According to lDrainel d2003l) . silicates of the type Mg-Fe py- 
roxenes (Mg x Fei_ x SiC<3) and olivines (Mg2 X Fe2-2xSiC<4) con- 
tribute a substantial fraction of the total mass of interstellar dust. 
These molecules are the main cause for the depletion of Si. In ad- 
dition, Mg, Fe, and O could be in oxide form (MgO, FeO, Fe2C<3, 



Silicates 


Composition 


A 


(x)MgSi0 3 , (l-x)FeSiO, 


B 


MgFeSi0 4 


C 


(x)Mg 2 Si0 4 , (l-x)Fe 2 Si0 4 


Others 


A+B, A+C, B+C, A+B+C 


Mg oxides 


Composition 


MgO 


Fe oxides 


Composition 


(x)FeO, (y)Fe 2 3 , (l-x-y)Fe 3 4 



Fe304). To estimate the amount of oxygen depletion in the Orion 
nebula, we investigate all combinations of dust grain composi- 
tions presented in Table [6] We implemented and IDL procedure 
that calculates Od us t for each of the combinations. The input pa- 
rameters of the program are the stellar and nebular Si, Fe, and 
Mg abundances, plus the nebular gas phase oxygen abundance. 
Once the total amount of Si, Fe, and Mg available to form dust 
grains is computed, the program determines the amount of O, Fe, 
and Mg used to form silicates (the whole amount of Si depleted 
is used to this aim); then, with the remaining Fe and Mg avail- 
able, the amount of O used to form Fe and Mg oxides is derived. 
Finally, the total amount of oxygen in the nebula (gas + dust) is 
determined for each combination of dust grain composition by 
adding the amount of oxygen depleted in silicates and oxides to 
the nebular gas phase abundance. 

We ran the IDL program 1000 times with random indepen- 
dent values for each of the input abundance^ to derive mean 
values and uncertainties of the amount of oxygen trapped in dust 
for each of the dust grain compositions considered. The results 
are presented in the upper panel of Fig. [2] We find that the oxy- 



5 Each of the input abundances has a normal distribution with mean 
values and standard deviations as those indicated in Tables [Tjand[5] 
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gen depletion is in the range 125-135 ppM, with corresponding 
uncertainties of ~ 10 pp M. This agrees very well wi th the value 
of 0.12 dex estimated bv lMesa-Delgado et all d2009l) . 

The total amount of oxygen in the nebula (gas + dust) de- 
termined for each combination of dust grain composition is pre- 
sented in the bottom panel of Fig. [2] Two cases are shown: the 
green boxes represent the results using CEL abundances, while 
the red boxes represent the results using RL abundances. We see 
that the correction for O depletion owing to dust grains does not 
depend much on the details of the dust composition in the dif- 
ferent combinations of dust grains we considerecQ. For compar- 
ison with the total nebular value, the stellar oxygen abundance 
is given in cyan. The height of the green and red boxes repre- 
sent the error bars combining the uncertainties in the gas phase 
abundance determination and in the dust phase abundance. The 
height of the cyan box represents the dispersion in the oxygen 
stellar abundances given in Table Q] The verdict from Fig. [2] is 
clear. Oxygen abundances in the Orion nebula derived from RLs, 
when corrected for the effects of depletion in dust grains, are 
very similar to those found in the B stars of the Ori OB 1 associa- 
tion. This is not the case for the oxygen CEL abundances, which 
lead to much lower values. We note that the oxygen CEL abun- 
dances would actually be compatible with the stellar abundances 
if using the intrinsic uncertainties of the analysis of individual 
stars (see Table[T]i. But as argued in Sect. [2] the degree of coher- 
ence between the abundances obtained for the 13 studied stars 
indicates that the oxygen abundance of the Ori OBI cluster is 
determined with much better accuracy than ~ 0.10 dex. 

Concerning the other elements considered in this study, there 
are observational arguments that a significant fraction of car- 
bon atoms must be locked in carbon-based dust and in PAHs 
although the a mount of carbon depletion in nebulae is still very 
poorly known dDraind l200l Uenkinsl 120091) . To our knowledge, 
there is no indication tha t nitrogen could be prese nt in solid form, 
except perhaps in ices dHilv-Blant et al.L 1201 Ot) . which are not 
expected to be present in ionized gas. As for neon, which is a 
noble gas, it is not expected to be incorporated in dust grains, at 
least not in significant amounts. 

5. Discussion 

Our study therefore leads us to the conclusion that the oxygen 
abundances derived in the Orion nebula from CELs are incom- 
patible with the abundances derived in the B stars (once the ef- 
fects of depletion in dust grains are taken into account), while 
RL oxygen abundances agree better (although the mean values 
seem to be somewhat higher that the stellar one). 

That RL abundances are more reliable than CEL abundances 
has been argued in many papers by Peimbert and coworkers (see 
e.g. iPeimbert & Peimbertl 120091) . including E04. One of the ar- 
guments is that the abundance discrepancy factor between RL 
and CEL abundances can be explained by temperature fluctua- 
tions of the same order as the difference between the temperature 
derived from the [O m] 4363/5007 line ratio and the one derived 
from the Balmer jump. In addition, temperature fluctuations of 
that order in the Orion nebula have been in ferred from the spa- 
tial variation of rj Oml across the nebula (lO'Dell et all 120031: 
iRubin et al.Lll993h . 



6 For completeness, we mention that Jenkins ( 2009) inferred from a 
detailed study of absorption line measurements in the solar vicinity that 
a significant proportion of interstellar oxygen should be trapped in ices. 
This statement, however, does not affect our estimates, since ices are 
not expected to survive in an ionized n ebula, because of the effect of 
photosputtering (Griaorieva et al., 2003). 
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Fig. 2 Upper panel: Amount of oxygen expected to be in solid 
form for various dust composition scenarios (see Table|6]). Lower 
panel: Comparison of the stellar oxygen abundance (wide cyan 
rectangle) with the gas+dust oxygen abundance in the nebula 
for the different dust scenarios. Green rectangles are obtained 
with CEL gas abundances, red rectangles are obtained with RL 
gas abundances. See text for an explanation of the considered 
uncertainties. 



There are, however, a certain number of controversial points 
to consider. 

First, when looking at Fig.[T] we see that the carbon RL abun- 
dance is slightly higher than the carbon abundance in the B stars, 
while we applied no correction for the depletion of carbon in 
dust grains. If we adopt the Si, Mg, and Fe depletions inferred 
from this figure, we infer a total dust-to-gas mass ratio in the 
silicates and oxides of 0.41+0.03%. With the "Orion mixture" 
considered in Cloudy, this implies a carbon depletion of ~50%. 
Then the total gas+dust carbon abundance, as derived from RLs, 
would be higher than the stellar one by about 0.2 dex. Note 
that the CEL carbon abundance is also higher than in the stars. 
However, in that case, the error bar is large and the CEL value 
is not really reliable. The nitrogen CEL abundance is slightly 
above the stellar one. However, taking into account the large icf 
(see Table|4| and the larger error bar on stellar N abundances (see 
Table Q]), both values are perhaps compatible. For neon, the RL 
abundance is somewhat uncertain and does not provide a useful 
constraint, but the CEL abundance is compatible with the stellar 
value if taking the model icfs and lower by about 0.2 dex if using 
the classical icfs. If the oxygen RL/CEL abundance discrepancy 
is caused by temperature fluctuations, one would expect the CEL 
abundances for C, N, and Ne to be significantly below the stel- 
lar values, and the RL abundances equal to the stellar ones. This 
is clearly not what occurs for carbon (for neon, the diagnostic 
depends on what icf is used). 

The second problem is related to the energy budget. To show 
this, we constructed a photoionization model with exactly the 
same parameters as the model shown in Fig. IB. 21 except that 
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for C, O, and Ne we take the RL abundances listed in Table [3] 
(and for N we also assume an RL-like abundance, by adding 
0.14dex to the CEL value). This model is shown in Fig. IB. 31 As 
seen in that figure, the value of r e [Om] returned by the model 
is below the observed one by 2 sigmas. In order for the model 
to reproduce the observed r e [Om] , one would need to consider 
extra heating, at least in some zones. Then, the computed inten- 
sities of the forbidden lines would be signficantly higher than ob- 
served. A model with slightly lower input abundances, in which 
the chemical composition of the nebular gas is exactly the same 
as found for the stars (but taking depletion into account) gives 
basically the same result: the nebular model becomes too cool 
and additional heating would make the mismatch between pre- 
dicted and observed forbidden line intensities even worse. 

Could it be, then, tha t the solution lies in o xygen-rich con- 
densations, as studied by IStasihska et al.l d2007l) ? This does not 
seem to be the case, at least in Orion. As found by those authors, 
oxygen rich condensations would bias the RL abundances up- 
wards , while we found in this paper that the RL abundances are 
compatible with the stellar ones. 

6. Conclusion 

The aim of the present paper was to re-examine the RL-CEL 
nebular abundance discrepancy problem in the light of recent 
high-quality abundance determinations in young stars. In this 
study, we focused on one object: the Orion star-forming region. 
There is no H n region in the solar neigborhood in which the stel- 
lar and gas phase abundances can be derived with better accu- 
racy, and the stars from the Orion OB 1 association, which share 
the same composition as the associated H n region, have been 
studied in great detail. 

We have re-evaluated the CEL and RL abundances of sev- 
eral elements in the Orion nebula and estimated the associated 
uncertainties, taking into account the uncertainties in the ioniza- 
tion correction factors for unseen ions. We estimated the amount 
of oxygen trapped in dust grains for several scenarios for dust 
formation. We compared the resulting gas+dust nebular abun- 
dances with the stellar abundances of a sample of 13 B-type 
stars from the Orion star-forming region (Ori OBI), analyzed in 
Papers I and III of this series. 

We find that the oxygen nebular abundance based on recom- 
bination lines agrees much better with the stellar abundances 
than that derived from the collisionally excited lines. 

However, there remain problems to be solved. First, the 
abundances of the other elements, C, N, and Ne, although ad- 
mittedly less accurate, do not deliver such a clear-cut picture. 
If the RL-CEL abundance dis crepancy were caused by tem pera- 
ture fluctuations, as argued bv lPeimbert & Peimbertl (12009). one 
should observe the same kind of bias in the CEL abundances of 
the other elements. Hopefully, the consideration of further ele- 
ments, such as S and Ar, whose abundances will soon be avail- 
able for the stars of the Orion OB 1 association, will bring use- 
ful additional constraints. Another problem is that with the RL 
abundances, the energy balance of the Orion nebula is not well 
understood. Investigating this problem will require detailed pho- 
toionization modeling and an accurate account of all heating and 
cooling processes taking place in this specific case, as well as 
the effect of uncertainties in the collision strenghts for the tran- 
sitions of interest. In parallel to such studies, another avenue to 
start exploring is the possibility of remaining biases in the stellar 
abundance determinations, for example because of a slight offset 
in the description of the ionization of oxygen in the model atmo- 
spheres of B-type stars, which unfortunately cannot be checked 



dire ctly here, since onl y On lines are seen in the studied stars 
rsee lSim6n-Dfazll2010l) . 

An important step forward would be to repeat the analysis 
performed in this paper for objects where the RL-CEL discrep- 
ancy is larger. The choice of the Orion nebula, while imposed 
by the quality of relevant observations and data analysis, is in- 
deed not the best one to investigate the RL-CEL nebular abun- 
dance discrepancy problem in H n regions because the discrep- 
ancy there is actually quite modest: ^ 0. 12 dex. With an adequate 
observing strategy and data handling, it should be possible to in- 
vestigate the RL-CEL abundance discrepancy in different envi- 
ronments, for example at low metallicities. This is a necessary 
condition to get a definite handle on this exasperating problem. 

Acknowledgements. Financial support from the Spanish Ministerio de Ciencia 
e Innovation under the project AYA2008-06166-C03-01 is acknowledged. 
This work has also been partially funded by the Spanish MICINN under the 
Consolider-Ingenio 2010 Program grant CSD2006-00070: First Science with the 
GTC I http://www.iac.es/consolider-ingenio-gtc I. 



References 

Aggarwal, K. M., & Keenan, F. P. 1999, ApJS, 123, 311 

Baldwin, J. A., Ferland, G. J., Martin, P. G., Corbin, M. R., Cota, S. A., Peterson, 

B. M., & Slettebak, A. 1991, ApJ, 374, 580 
Blagrave, K. 2006, Ph.D. Thesis 

Bresolin F, Gieren W., Kudritzki R.-R, Pietrzyhski G., Urbaneja M. A., Carrara 

G., 2009, ApJ, 700, 309 
Bresolin F, Stasinska G, Vflchez J. M., Simon J. D., Rosolowsky E., 2010, 

MNRAS, 404, 1679 
Butler, K., & Giddings, J. R. 1985, in Newsletter of Analysis of Astronomical 

Spectra, No. 9 (Univ. London) 
Charbonnel, C, Schaerer, D., & Meynet, G. 2003, Astronomical Society of the 

Pacific Conference Series, 304, 
Cunha, K., & Lambert, D. L. 1992, ApJ, 399, 586 
Cunha, K., & Lambert, D. L. 1994, ApJ, 426, 170 
Draine, B. T. 2003, ARA&A, 41, 241 
Ercolano B., 2009, MNRAS, 397, L69 

Esteban, C, Peimbert, M., Torres-Peimbert, S., & Escalante, V. 1998, MNRAS, 
295, 401 

Esteban, C, Peimbert, M., Garcfa-Rojas, J., Ruiz, M. T., Peimbert, A., & 

Rodriguez, M. 2004, MNRAS, 355, 229 
Ferland, G. J., Korista, K. T, Verner, D. A., Ferguson, J. W., Kingdon, J. B., & 

Verner, E. M. 1998, P ASP, 110, 761 

Froese Fischer 2004, http://hf8.vuse.vanderbilt.edu 

Froese Fischer C, Tachiev G., 2004, ADNDT, 87, 1 

Galavis, M., Mendoza, C, Zeippen, C. J., 1995, A&AS 111, 347 

Galavis, M.E., Mendoza, C, Zeippen, C.J., 1997, A&AS 123, 159 

Galavis, M. E„ Mendoza, C, Zeippen, C. J., 1998a, A&AS 131, 499 

Galavis, M. E., Mendoza, C, Zeippen, C. J., 1998b, A&AS 133, 245 

Garcfa-Rojas, J., & Esteban, C. 2007, ApJ, 670, 457 

Giddings, J. R. 1981, Ph.D. Thesis, University of London 

Grigorieva, A., Thebault, P., Artymowicz, P., & Brandeker, A. 2007, A&A, 475, 

755 

Hillier, D. J., & Miller, D. L. 1998, ApJ, 496, 407 

Hily-Blant, P., Walmsley, M., Pineau Des Forets, G., & Flower, D. 2010, A&A, 
513, A41 

Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 

Hudson, C. E., & Bell, K. L., 2004, MNRAS 348, 1275 

Jenkins, E. B. 2009, ApJ, 700, 1299 

Keenan, F. P., Aller, L. H., Bell, K. L., Hyung, S., McKenna, F. C, & 

Ramsbottom, C. A. 1996, MNRAS, 281, 1073 
Kurucz, R. L. 1993, CD-ROM No. 13 (Cambridge, Mass.: SAO) 
McLaughlin B. M., Bell K. L., 2000, JPhB, 33, 597 
Mendoza, C, Zeippen, C. J., 1982a, MNRAS 202, 981 
Mendoza, C, Zeippen, C. J., 1982b, MNRAS 198, 127 
Mendoza, C, Zeippen, C. J., Storey, P. J., 1999, A&AS 135, 159 
Mendoza, C, 1983, IAUS, 103, 143 

Menten, K. M., Reid, M. J., Forbrich, J., Brunthaler, A. 2007, A&A, 474, 515 
Mesa-Delgado, A., Esteban, C, & Garcfa-Rojas, J. 2008, ApJ, 675, 389 
Mesa-Delgado, A., Esteban, C, Garcia-Rojas, J., Luridiana, V., Bautista, M., 

Rodriguez, M., Lopez-Martin, L., & Peimbert, M. 2009, MNRAS, 395, 855 
Astrophysics of gaseous nebulae and active galactic nuclei, 2nd. ed. by D.E. 

Osterbrock and G.J. Ferland. Sausalito, CA: University Science Books, 2006 



S. Simon-Diaz & G. Stasinska: Stars, gas and dust in Ori OB 1 



Table 7 Atomic data for line emission used in Sect. [3] 



collision strengths 



N II Hudson & Bell (2004) + priv. 



1999) 
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Ne III McLaughlin & Bell (2000 



SII 
S III 
Ar III 
ArlV 



transition probabilities 
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Taval & Gupta (1999) 

Galavis, Mendoza & Zeippen ( 1998b) 

Ramsbottom. Bell & Keenan (1997) + priv. 
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Storey & Zeippen (2000) 
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Mendoza. Zeippen & Storey (1999) 
Storey & Zeippen (2000) 
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Appendix A: Atomic data for nebular emission line 
Appendix B: Photoionization model 

To estimate the ionization correction factors as well as possi- 
ble for the elements whose ionic abundances were derived from 
the optical spectrum of E04, we constructed a photoionization 
model reproducing the observed characteristics of this spectrum. 

The model was constru cted with the photoi onization code 
Cloudy (as last described in lFerland et allTl998l) . using version 
8.0. We considered an open geometry (as appropriate for the 
Orion nebula) and taken the "intrinsic line fluxes" and "aperture 
slit" options. Of the options offered by Cloudy, this one seems 
the most approriate for the E04 observation. 

Because we are here mainly interested in the icfs, we 
do not try to fit all observed line intensities (we already 
know that this would not be possible for RLs and CELs 
simultaneously). On the other hand, we aim at reproduc- 
ing the line ratios indicative of density ([Ariv] 4740/4711, 
[On] 3726/3729, and [Sn] 6731/6716) and ionization struc- 
ture ([Oiii] 5007/[On] 3727, [Sm] 9069/([Sn] 6731+16) and 
[Ariv] 4711+40/[Arm] 7135). 

We assumed an exponential density la w, which is usuall y 
considered for the Orion nebula (see e.g. iRubin et"ail Il993h : 



N(H) - N ma _ x exp[(x-x p )/h], with N(H)=N nmx for x > x p . We 
considered as free parameters h, N mix , and x p , and fixed the in- 
ternal radius of the model to be ro = 0.9 x p . 

The second important input to the photoionization model is 
the ionizing radiation field. We considered the stellar parame- 
ters derived for 6 1 Ori C by ISimon-Dfaz et all d2006l) and ob- 
tained the spectral energy distribution (SEP) as predicted b y 
the stellar atmopshere codes WMb asic dPauldrach et all 12 001). 
CMF GEN dHillier & Milieilll998l). T LUSTY dHubenv & Lanzl 
1 1995b . and FASTWIND dPuls et ail 120051) . Fig. ED compares 
the four SEDs, also indicating the ionization energies for vari- 
ous of the ions present in the O rion nebula. As commented in 
ISimon-Dfaz & Stasinska! d2008h . there are important discrepan- 
cies between the ionizing SEDs predicted by the four stellar at- 
mosphere codes. We considered the four SEDs in our photoion- 
ization models, and it turned out that WMbasic gives the best 
solution for this case (this argument is based on the simultane- 
ous fitting of the ionization degree constraints, corresponding to 
the ionization of S + , + , and Ar 2+ ; see also Figure lBTTT i. We also 
explored slight changes in the effective temperature, and found 
that the indicated stellar parameters (T e ff = 39000 K, log g_ =4.1) 
provide the best solution. 

In all cases, a total number of hydrogen ionizing photons, 
Q(H°)=10 48 67 photons s was considered. This value was ob- 
tained assuming the sam e m v magnitude and extinction as in 
ISimon-Dfaz et alJ d2006l) . but adopting a distance to the Orion 
nebula of 400 pqj, instead of the value originally adopted 450 
pc. 

Finally, we took grains into account in the computation of 
all the models, assuming the "Orion type" grain composition 
considered by Cloudy, scaled by a factor 1.1. This scale factor 
was included to obtain the same abundance of Si, Mg, and Fe in 
grains as computed in this paper (see Sect. |4j. 

To judge the quality o f the model, we use the same procedure 
as IStasiriska et alJ d2010l) . plotting for each observable O (i.e. 
here for each relevant line ratio) the value of 

k(0) = (logO mod - logO obs )/r«9), (B. 1) 

where <9 mo d is the value returned by the model, <9 b s is the ob- 
served value, and t(O) the accepted tolerance in dex for this ob- 
servable. For each observable, the value of t(O) is defined by 

r(0) = log(l+er(0)), (B.2) 



7 This value agrees b etter with recent dete rminations of the distance 
to the Orion nebula (seeM enten et all l2007. and references therein) 
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Fig. B.l Comparison of the ionizing SEDs for 1 OriC as pre- 
dicted by various stellar atmospher e models. Stellar paramet ers: 
T eff =39000, logg_=4.0 (derived by ISimon-Diaz et all [20061 by 
means of a spectroscopic analysis of the optical spectra using 
FASTWIND) 



modelQ We may notice that using the CEL abundances implies a 
T e [0 m] in the model much higher than observed. This suggests 
that the input abundances (especially that of oxygen) might be 
too small. Since we argue in Sect. [5] that the gas phase oxygen 
abundance seems to be better estimated by RLs than by CELs, 
one may wonder whether it is correct to fit the ionization state 
of the nebula with ratios of CELs, as we have done. We believe 
that this procedure is reasonably correct, unless the discrepancy 
between RL and CEL abundances is much different in the low 
ionization and in the high ionization zone. 

Figure IB. 31 shows the results of another model where this 
time the input abundances of C, O, and Ne are those obtained 
from the RL lines and the abundance of N is the CEL one in- 
creased by 0.14 dex. The aim of this model is to discuss the 
energy budget problem as presented in Sect. [5] 



Table B.l List of nebular line ratios used to compare with the output 
of the photoionization models and the corresponding observed values 
(from E04). 



Parameters 


Line ratio 


Oobs 


e,(0) 


Density (constraint) 


A^Ariv] 


[Ariv] 4740/4711 


1.21 


0.08 


We [On] 


[On] 3726/3729 


2.07 


0.06 


2V e [Sn] 


[Sn] 6731/6716 


1.81 


0.08 


Ionization degree (constraint) 


s i+ /s + 


[S in] 9069/[Sn] 6716+31 


5.53 


0.16 


o 2+ /o + 


[Om] 5007/[On] 3727 


4.64 


0.06 


Ar 3+ /Ar 2+ 


[ Ar rv] 47 1 1 +40/[ Ar m] 7135 


0.0136 


0.11 


Abundances 


O 


100X [Om] 5007/Hg 


383 


0.03 




100X Oii4075/Hs 


0.27 


0.10 


C 


lOOx Cn4267/rL; 


0.23 


0.05 


N 


lOOx [Nn] 6584/FL; 


37.7 


0.05 


S 


100X [S m] 9069/H/j 


30.2 


0.14 


Ar 


100X [Arm] 7135/H^ 


16.2 


0.07 


Ne 


100X [Nem] 3869/Hs 


22.9 


0.04 


Fe 


100X [Fern] 4659/fy 


0.55 


0.05 


Temperature 


rjOm] 


[Om] 4363/5007 


0.00339 


0.03 


r e [Sm] 


[S m] 6312/9069 


0.0613 


0.15 


r e [On] 


[On] 7325/3727 


0.116 


0.13 


r e [Nn] 


[Nn] 5755/6584 


0.0180 


0.06 



where e r (0) is the maximum "acceptable" relative error on the 
observable. A model is satisfactory only if for each of the ob- 
servables used as constraints, k(0) is found to be between -1 
and + 1 . Table lB.ll lists the values of the observables used as con- 
straints to fit the model and the corresponding values of e r (0) at- 
tributed to them (it also lists, for information, some other observ- 
ables that are not used as constraints, such as line ratios indica- 
tive of the electron temperature or of the element abundances). 

Figure IB. 21 is a graphic representation of one of our best 
models and of the quality of the fits. In this model (Model B.l) 
the input abundances are the CEL ones given in Table [5] In the 
right panel of the figure the icfs are listed corresponding to the 



8 In Tables [4] and [3] are listed the values of icfs that were used to 
compute the element abundances. We estimated the uncertainties in the 
icfs by changing the input parameters of the photoionization model, still 
keeping the values of k in the acceptable range. 
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Fig. B.2 Results from a photoionization model in which the input abundances are the CEL ones given in Table[3] The numbers to the 
right of each line ratio indicate the deviation from the observed value in units of k{0) and the value obtained from the photoionization 
model (the observed value being given in Table IB. ll . In the right column are listed the logarithms of the icfs corresponding to this 
model. 
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Fig. B.3 Results from a photoionization model in which the input abundances are the RL ones given in Table[5] The numbers to the 
right of each line ratio indicate the deviation from the observed value in units of k(0) and the value obtained from the photoionization 
model (the observed value being given in Table IB. ll . In the right column are listed the logarithms of the icfs corresponding to this 
model. 



